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1.  Introduction 


As  the  electronics  behind  the  future  warrior  systems  become  more  sophisticated,  the  weight  of 
the  batteries  is  an  ever-increasing  burden.  One  solution  is  to  create  a  compact  micro-burner 
device  as  shown  in  figure  1  (Norton  et  ah,  2004)  that  combusts  a  higher  energy  density  fiiel  such 
as  methane  (energy  density  =  3053  W-hr/kg  compared  to  125  W-hr/kg  for  lithium  (Li)-ion 
batteries)  and  converts  the  released  enthalpy  into  electrical  power. 


Figure  1.  Micro-burner  schematic. 

Note:  The  catalytic  insert  consists  of  nanosized 

platinum  (Pt)  dispersed  within  porous  alumina  (AI2O3). 

Although  there  are  many  computational  studies  that  detail  the  complete  hydrogen  (H2) 
combustion  mechanism  for  reactant  and  product  species  interacting  with  the  catalytically  active 
Pt  cluster  (Mhadeshware  and  Vlachos,  2007),  few  studies  consider  the  effect  of  the  AI2O3 
support.  New  reaction  pathways  can  arise  due  to  support  surface  termination  and  reactivity  at 
the  Al203/Pt  interface.  One  example  of  such  a  pathway  is  the  “inverse  spillover  effect”  (ISE), 
which  occurs  when  water  (H2O)  chemisorbs  or  dissociates  on  the  support  forming  mobile  species 
that  can  migrate  to  the  catalytically  active  particle  and  further  promote  combustion. 

Experimental  evidence  for  ISE  comes  from  the  work  of  Wang  et  al.  (1996)  who  demonstrated 
that  carbon  monoxide  (CO)  can  liberate  H2  from  an  H2O  bound  to  AI2O3  support.  In  this  report, 
we  seek  to  augment  the  current  model  for  H2  micro-combustion  to  include  reactions  that  are 
initiated  on  the  catalytic  support.  Specifically,  we  propose  a  mechanism  by  which  H2O  and  H2 
dissociate  on  the  support  forming  hydroxyl  groups,  which  can  further  dissociate  and  diffuse  to 
the  catalytic  particle. 
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2.  Model  Parameters  and  Validation 


2.1  Alpha-alumina  oxide  (aA^Os)  Surface  Termination 


Our  model  (figure  2)  consisted  of  a  nine  atomic  layer  thick,  aluminum  (Al)-terminated  (0001) 
slab  that  is  repeated  under  periodic  boundary  conditions  as  a  2x2  supercell  with  PI  symmetry 
and  a  30A  vacuum  layer  to  prevent  any  interaction  between  periodic  images.  We  chose  this 
surface  based  on  the  availability  of  experimental  and  theoretical  data  in  the  literature  as  well  as 
the  work  of  Marmier  et  al.  (2004),  who  calculated  surface  phase  diagrams  as  a  function  of 
temperature  and  the  oxygen  (O2)  and  H2  partial  pressures  for  several  different  crystal  orientations 
and  surface  terminations.  The  lattice  parameters  for  the  rhombohedral  unit  cell  (a=b=  4.749  A, 
c=12.991  A)  were  taken  directly  from  experimental  results  (Swansen,  1960)  and  were  not 
optimized  during  the  simulation.  In  addition,  we  constrained  the  bottom  two  layers  of  the  slab  to 
reflect  the  bulk  AI2O3  geometry. 


Figure  2.  Cross  section  of  fully  relaxed  Al-terminated 
aAl203  (0001)  slab  used  for  calculation. 

The  O  and  Al  atoms  are  red  and 
magenta,  respectively. 


2.2  Model  Parameters 

All  calculations  were  performed  using  ideal  conditions  (0  K,  ultra  high  vacuum,  defect  free 
surface).  The  calculations  were  executed  within  the  DMol  (Delley,  2000)  module  of  the 
Materials  Studio  (version  4)  software  package  using  a  double-numeric  basis  set  with  polarization 
functions  (DNP)  and  the  Perdew-Burke-Emzerhof  (PBE)  (Perdue,  1996)  version  of  the 
generalized  gradient  approximation  (GGA)  to  represent  the  electron  exchange  and  correlations. 
The  ion  cores  were  described  by  a  density  functional  semi-core  pseudopotential  (DSPP)  (Delley, 
2002). 

To  validate  our  calculations,  we  compare  our  results  for  the  surface  reconstruction  of  the  relaxed 
aAl203  (0001)  slab  with  the  results  of  other  theoretical  and  experimental  investigations  (table  1). 
In  agreement  with  other  theoretical  studies,  our  simulation  predicts  an  89%  contraction  of  the 
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inter-atomic  spacing  of  top  surface  layer  and  6%  expansion  of  the  first  sub-layer  for  the  ultra¬ 
clean  aAl203  (0001)  surface.  The  predicted  surface  reconstruction  was  explained  by  Sousa  et  al 
(1993)  as  being  a  result  of  charge  redistribution  due  to  the  highly  ionic  nature  of  alumina.  The 
experimental  value  for  this  relaxation  is  closer  to  ~50%.  The  discrepancy  between  theoretical 
prediction  and  experimental  measurements  may  be  due  to  the  difficulty  in  preparing  a  perfectly 
terminated  surface  with  no  adsorbed  atoms  or  defects. 

Table  1.  Comparison  of  theoretical  and  experimentally  measured  changes  in  the  inter-atomic  layer  spacing  of  Al- 
terminated  aAhOs  (0001)  slab  with  respect  to  their  unrelaxed  geometry. 


Theoretical 

Experimental 

Ours 

Hinneman 

Verdozzi 

Hass 

Alavi 

Ruberto 

Carrasco 

Guenard 

Ahn 

Wxygen 

layers 

3 

9 

18 

3 

3 

9 

11 

Functional 

PBE/DSPP 

PBE/PAW 

LDA/NCPP 

PBE/NCPP 

PW91/USPP 

PW91/NCPP 

PBE/PAW 

-89.2 

-86.4 

-87.4 

-98 

-97 

-85.5 

-93.8 

-51 

63 

+6 

+4 

+3.1 

+5 

+2 

+3.2 

+6.1 

+16 

Ail3J.Af^J 

-39.9 

-45.4 

-41.7 

-48 

-53 

-45.4 

-46.1 

-29 

A114J.OL5J 

+18.9 

+20.5 

+18.3 

+21 

-27 

+19.8 

+22.0 

+20 

QlSJ.Af^J 

+17.1 

+5 

+5.6 

+4.8 

+8.5 

-31.2 

-6.8 

-8.3 

-7.1 

-11.6 

Ailvj.oL«J 

0 

+1.3 

+1.1 

+1.3 

+2.2 

0 

-1.3 

-0.5 

-0.8 

+0.7 

Aim.Af^J 

+4.6 

+6.4 

+3.0 

+3.8 

Ail'»J.o^"J 

-1.2 

-0.6 

-0.7 

-3.2 

To  simulate  adsorption  phenomena,  we  added  one  adsorbate  molecule  per  supercell,  which  is 
equivalent  to  approximately  1/12  monolayer  according  to  Verdozzi  et  al.  (1999),  who  define  a 
monolayer  as  having  one  metal  atom  per  surface  oxygen.  Binding  energies  were  computed  by 
subtracting  the  energy  of  the  clean  fully  relaxed  slab  and  the  adsorbent  molecule  (H2O,  O2,  H2) 
from  the  total  energy  of  the  system  after  adsorption.  We  performed  barrier  calculations  using  the 
linear  synchronous  transit  (LST)  method  of  Govind  et  al.  (2003)  to  extrapolate  between  the 
reactant  and  product  structures  along  the  diffusion  pathways  illustrated  in  figure  3. 
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Figure  3.  Top  view  of  a  2x2  supereell  showing  allowed  surfaee 

binding  sites.  The  numbers  correlate  with  the  dissociation 
product  and  diffusion  path  notation  referred  to  throughout 
the  report.  When  referring  to  dissociation  products,  the 
endpoints  of  the  line  indicate  the  atoms  to  which  the 
dissociated  species  bind. 


2.3  Electron  Spin  State 

The  reaction  chemistry  of  oxygen  involved  reactions  on  AI2O3  surfaces  cannot  be  adequately 
described  without  careful  consideration  of  the  triplet-to-singlet  spin  conversion  that  occurs  when 
the  IpTIg*  orbitals  hybridrize  with  the  surface  states.  We  find  the  energy  difference  between  the 
triplet  and  the  singlet  ground  state  of  free  O2  to  be  19.3  kcal/mol,  which  is  in  good  agreement 
with  the  experimentally  measured  value  of  22.6  kcal/mol  (Herzberg,  1950).  Figure  4  shows 
adsorption  structures  and  dissociation  products  for  O,  H,  O2,  H2,  and  H2O  on  the  Al-terminated 
aAbOs  (0001)  surface.  For  all  adsorption/dissociation  products  except  figure  4A,  4E,  and  4J,  the 
lowest  energy  spin  state  was  singlet.  In  table  2,  we  detail  the  key  molecular  features  for  these 
structures  as  well  as  their  binding  energies.  From  table  2,  it  is  clear  that  triplet  states  result  in 
tighter  0=0  bonds  and  a  local  elongation  of  the  Al'^'^-O  bonding  scheme.  The  spin  state  has  no 
apparent  effect  on  the  length  of  hydroxyl  or  bonds.  Based  on  the  binding  distances,  it 

would  appear  that  the  triplet  to  singlet  spin  conversion  occurs  between  2  and  1.9A  above  the 
surface  for  molecular  O2  and  between  1.8  and  1.5  A  for  atomic  oxygen.  However,  this  can  only 
be  verified  by  carefully  sampling  both  the  singlet  and  triplet  potential  energy  surfaces  as  the 
adsorbing  species  approaches  the  surface. 
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Figure  4.  Adsorption  and  dissociated  structures  for  (A)  oxygen  tetrahedron,  (B)  oxygen  bridge,  (C)  hydrogen 
tetrahedron,  (D)  hydrogen,  (E)  1-1  molecularly  adsorbed  O2,  (F)  1-2  dissociated  O2,  (G)  2-2 
dissociated  O2,  (H)  1-2  dissociated  H2,(I)  1-4  dissociated  H2,  (J)  2-2  dissociated  H2,(K)  1-1 
molecularly  adsorbed  H2O,  (L)  1-2  dissociated  H2O,  (M)  1-4  dissociated  H2O,  and  (N)  the  key 
indicating  atomic  layer  to  which  atom  originally  belonged. 
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Table  2.  Effect  of  spin  state  on  O2  adsorption  binding  energies  and 
bond  lengths.  The  ID  notation  refers  to  the  structures  in 
figure  4. 


Atomic  Oxygen  Adsorption 

^binding 

Q[ads]_Q[2] 

Q[ads]_^j[l] 

^bond 

(kcal/mol) 

(A) 

(A) 

(°) 

Triplet  A 

[-42]" 

1.782 

112 

Singlet  A 

-30  [-35] " 

1.767 

115 

Triplet  B 

r-isr 

1.508 

1.797 

49 

Singlet  B 

-49  [-53] " 

1.546 

1.803 

50 

Atomic  Hydrogen  Adsorption 

^binding 

jj[ads]_Q[2] 

jj[ads]_^j[l] 

^bond 

Triplet  C 

-20 

1.626 

114 

Singlet  C 

-36 

1.625 

114 

Triplet  D 

-98 

0.971 

124 

Singlet  D 

-116 

0.972 

124 

Molecular  Oxygen  Adsorption  /  Dissociation 

^binding 

0=0 

Q[ads]_^j[l] 

^bond 

Triplet  E 

-7 

1.255 

1.997 

107 

Singlet  E 

-13 

1.276 

1.959 

110 

Triplet  F 

21 

1.364 

1.895 

74 

Singlet  F 

-3 

1.393 

1.851 

76 

Triplet  G 

51 

1.499 

1.837 

48 

Singlet  G 

25 

1.522 

1.823 

49 

Molecular  Hydrogen  Dissociation 

^binding 

jj[ads]_Q[2] 

jj[ads]_^j[l] 

^bond 

OH,HOAl 

Triplet  H 

-15 

0.982 

1.594 

110,  120 

Singlet  H 

-14 

0.981 

1.593 

110,  120 

Triplet  I 

-14 

0.971 

1.621 

120,  113 

Singlet  I 

-10 

0.972 

1.615 

121,  112 

Triplet  J 

-85 

0.973 

127,  — 

Singlet  J 

-67 

0.973 

127,  — 

H2O  Adsorption  /  Dissociation 

^binding 

jj[ads]_Q[2] 

^OAIO 

OH,OA10 

Triplet  K 

-27  [-23]  ^ 

0.982 

1.982 

— ,  87 

Singlet  K 

-26  [-23]  ^ 

0.983 

1.987 

— ,  86 

Triplet  F 

-40  [-33]  ^ 

0.982 

1.740 

110,  112 

Singlet  F 

-38  [-33]  ^ 

0.981 

1.740 

109,  114 

Triplet  M 

-38  [-33]  ^ 

0.972 

1.759 

121,97 

Singlet  M 

-34  [-33]  ^ 

0.973 

1.758 

122,  97 

^The  number  in  brackets  is  from  Gamallo,  2007 
'’The  number  in  brackets  is  from  Haas,  2000. 


Since  an  objective  of  this  research  is  to  establish  whether  dissociated  oxygen  from  the  alumina 
support  can  diffuse  to  the  catalytic  particle  and  all  possible  surface  diffusion  pathways  include  a 
bridging  conformation,  we  chose  to  perform  our  barrier  calculations  using  a  singlet  spin  state  for 
both  the  reactant  and  products.  This  is  the  most  accurate  method  to  explore  the  effects  of  crystal 
symmetry  and  surface  reconstructions  on  the  dissociation  and  diffusion  barriers  for  reaction 
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pathways  that  do  not  involve  a  spin  change.  In  future  studies,  we  will  perform  a  detailed 
sampling  of  both  the  triplet  and  singlet  potential  energy  surfaces  using  the  constrained  geometry 
method  to  assess  the  reaction  barriers  for  pathways  that  involve  spin-to-triplet  transformations 


3.  Results  and  Discussion 


3.1  Effects  of  Adsorption  on  Surface  Reconstruction 

In  table  3,  we  demonstrate  the  effect  of  the  adsorption  of  different  species  on  both  the  local  and 
long-range  reconstruction  of  the  aAl203  (0001)  surface.  Both  H  and  O  can  directly  bind  with 
similar  binding  energies  (~  -30  kcal/mol)  to  the  surface  Al'^'^  forming  a  tetrahedron  with 
neighboring  atoms.  By  saturating  the  surface  Al'^*\  the  0“'^^  reduces  the  driving  force  for  the 
contraction  of  the  first  inter-atomic  layer.  Although  H  binds  closer  to  the  surface  Al'^'^  than  O, 
both  H  and  O  produce  the  same  elongation  of  the  local  Al'^'^-O'^^^  bonding  scheme.  H  and  O  can 
also  directly  bind  to  the  atoms.  However,  in  this  configuration,  there  are  significant 
differences  in  both  the  adsorption  structures  and  binding  energies  (Ebind^  ~  -50,  Ebind^  ~  - 
1 16kcal/mol).  Whereas  O  pulls  the  Al'^'^  ~  8%  away  from  the  surface,  H  drives  Al'^'^  deeper  into 
the  lattice.  The  closer  the  binding  site  is  to  an  Al'^'^  site,  the  deeper  it  is  driven  into  the  lattice  and 
the  greater  the  asymmetrical  lengthening  of  the  local  bonding  scheme.  This  has  a 

profound  effect  on  the  relaxation  of  the  top  two  inter-atomic  layers.  As  the  surface  Al'^'^  is  driven 
into  the  lattice,  it  reduces  the  ability  of  the  oxygen  layer  to  charge  compensate  for  excess  charge 
on  the  remaining  unsaturated  surface  aluminum  atoms,  resulting  in  a  smaller  contraction  of  the 
first  inter-atomic  layer. 
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Table  3.  Effect  of  adsorption  on  surface  reconstruction. 


Layerl 

Layer2 

AAl'*' 

Allil-o'^l 

(%) 

(%) 

(%) 

A 

A1203 

-89 

+6 

1.704 

1.704 

1.704 

0(A) 

-83 

+6 

+  11 

1.783 

1.783 

1.783 

0(B) 

-98 

+9 

+8 

1.826 

1.715 

1.714 

H(C) 

-87 

+6 

+11 

1.784 

1.784 

1.784 

H(D) 

Non  uniform 

-14 

1.841 

1.747 

1.746 

02(E) 

-90 

+7 

+7 

1.735 

1.735 

1.733 

02(F) 

-84 

+4 

+9 

1.878 

1.711 

1.710 

02(G) 

-86 

+13 

+13 

1.869 

1.862 

1.733 

H2(H) 

-83 

+10 

+13 

1.920 

1.764 

1.764 

H2(I) 

-84 

+9 

+11 

1.796 

1.790 

1.770 

H2(J) 

Non  uniform 

-16 

1.844 

1.837 

1.745 

H2O  (K) 

-104 

+6 

+6 

1.733 

1.721 

1.720 

H2O  (L) 

-83 

+9 

+13 

1.894 

1.755 

1.755 

H2O  (M) 

-85 

+7 

+11 

1.793 

1.779 

1.766 

Note:  Layer  1  and  Layer  2  are  the  percent  change  in  the  L*  and 
2"''  inter-atomic  layers  with  respect  to  their  bulk 
coordinates.  A  negative  sign  for  Layer  lor  2  represents  a 
contraction  of  the  layer.  calculates  how  much  the 

adsorbate  pulls  the  surface  Al^'^  site  from  its  original 
relaxed  position.  Al^'^-O^^'  are  the  surface  bonds 
neighboring  the  adsorption  site. 

O2  and  H2O  can  molecularly  adsorb  to  surface  with  binding  energies  of  ~20keal/mol; 
whereas,  H2  cannot.  O2  adsorbs  closer  to  the  surface  than  H2O  and  does  not  ehange  the 
contraction  of  the  first  inter-atomie  layer.  In  eontrast,  moleeularly  adsorbed  H2O  eauses  the 
surfaee  Al'^'^  to  eontraet  below  the  atoms,  changing  the  surface  termination  from  Al- 
terminated  to  0-terminated  although  the  A1  and  O  atoms  are  nearly  co-planar. 

As  shown  in  figure  4,  there  are  three  unique  eonfigurations  for  the  dissoeiated  products, 
henceforth  referred  to  as  1-2,  1-4,  and  2-2  dissoeiation.  H2  ean  form  all  three  dissoeiation 
produets;  however,  H2O  cannot  form  2-2  dissociation  products  and  O2  cannot  form  1-4 
dissociation  products.  In  comparing  the  dissoeiation  products,  the  following  trends  are  elear:  (1) 
dissoeiation  reduces  the  contraction  of  the  first  inter-atomie  layer  regardless  of  which  species  is 
dissoeiating  and  (2)  EbmdOi  >  EbindH2  >  EbindH20.  In  regards  to  whieh  type  of  dissociation 
product  has  the  lowest  energy,  it  depends  on  which  species  are  present.  For  both  H2O  and  O2, 
the  lowest  energy  dissociation  products  are  1-2;  however,  for  H2,  the  lowest  energy  dissociation 
product  is  2-2. 

3.1  O2  Molecular  and  Dissociative  Adsorption 

The  moleeular  and  1-2  dissoeiative  adsorption  of  O2  from  free  O2  appears  to  be  spontaneous. 
However,  the  actual  barrier  for  direct  and  indirect  (i.e.,  from  1-2  molecularly  adsorbed) 
dissoeiation  eannot  be  determined  using  the  LST  transition  state  method  because  the  reaction 
involves  a  triplet-singlet  spin  eonversion,  whieh  can  only  be  aceurately  described  by  sampling 
both  the  triplet  and  singlet  potential  energy  surfaees.  Nevertheless,  we  can  accurately  calculate 
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the  barriers  to  the  further  dissoeiation  of  1-2  into  2-2  produets  since  the  minimum  energy 
structure  for  both  of  these  structures  is  a  singlet  spin  state.  We  find  that  the  further  dissociation 
of  1-2  adsorption  products  is  highly  endothermic  and  not  kinetically  favorable  (Erxn  =  28 
kcal/mol,  Ebarrier=  53  kcal/mol).  Given  these  results,  the  indirect  dissociation  of  O2  is  not 
considered  a  viable  source  for  isolated  substrate  bound  atomic  oxygen. 

3.2  H2O  Adsorption  and  Dissociation 

The  question  arises  as  to  whether  dissociated  H2O  can  serve  as  a  source  of  substrate  bound 
oxygen.  To  answer  this  question,  we  investigated  the  following  five  pathways  for  the  adsorption 
and  dissociation  of  H2O: 

1.  Free  H2O  molecularly  (figure  4K)  adsorbs  to 

2.  Free  H2O  dissociates  into  a  1-2  conformation  (figure  4L)  by  forming  hydroxyls  with  Al'^’^ 
and  its  first  nearest  neighboring 

3.  Free  H2O  dissociates  into  a  1-4  conformation  (figure  4M)  by  forming  hydroxyls  with  Al'^’^ 
and  its  second  nearest  neighboring 

4.  Adsorbed  H2O  (figure  4K)  dissociates  into  a  1-2  conformation. 

5.  Adsorbed  H2O  (figure  4K)  dissociates  into  a  1-4  conformation. 

Our  results  indicate  that  H2O  spontaneously  dissociates  into  both  1-2  and  1-4  conformations. 
Both  dissociated  products  have  equivalent  adsorption  energies.  H2O  can  also  adsorb  molecularly 
to  surface  A1  with  a  slightly  higher  energy  than  its  dissociated  product  (~+10kcal/mol).  The 
barrier  for  the  molecular  adsorption  of  H2O  is  0.083  kcal/mol,  which  is  within  the  limits  of  the 
accuracy  of  the  model.  Once  molecularly  adsorbed,  the  barrier  to  further  dissociation  into  1-2 
and  1-4  products  are  23  and  9.6  kcal/mol,  respectively.  Experimentally  (Elam,  1998),  there  has 
been  no  evidence  found  for  the  existence  of  molecularly  adsorbed  H2O.  Since  the  energy 
released  due  to  spontaneous  dissociation  from  free  H2O  is  well  in  excess  of  the  barriers  for 
dissociation  from  molecular  H2O,  our  results  suggest  that  if  molecular  H2O  exists  on  the  surface, 
it  is  unlikely  to  have  a  long  lifetime. 

The  next  series  of  questions  that  arise  is  (1)  Can  hydrogen  diffuse  away  from  dissociated  H2O? 
and  (2)  Can  the  remaining  Al'^'^OH  hydroxyl  group  further  dissociate  into  a  H  and  O  pair?  We 
find  that  the  easiest  paths  for  the  H  to  diffuse  away  from  dissociated  H2O  are  paths  10,  13,  15 
(Ebarrier  ~  17-24  kcal/mol).  These  paths  avoid  the  influence  of  Al'^'^  and  atoms.  For  all  other 
paths  considered,  the  barrier  to  the  diffusion  of  H  from  dissociated  H2O  was  >40  kcal/mol.  The 
barrier  to  the  further  dissociation  of  the  remaining  Al'^’^OH  hydroxyl  into  2-2  and  2-4  0*»H 
dissociation  pairs  is  45  and  34  kcal/mol,  respectively. 
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3.3  O  Surface  Diffusion 


The  final  question  to  be  answered  is  whether  the  isolated  substrate  bound  O  that  is  produced 
from  the  dissociation  of  H2O  can  diffuse  across  the  AI2O3  surface  to  the  catalytic  particle  to 
promote  combustion.  In  table  4,  we  provide  barrier  calculations  for  diffusion  of  O  and  H 
between  all  allowed  surface  sites.  We  find  that  O  can  easily  diffuse  from  Al'^*^  tetrahedron 
(figure  4A)  to  a  bridging  (figure  4B)  conformation  (Ebarrier  =  4  kcal/mol).  Once  in  a  bridging 
conformation,  O  can  diffuse  between  bridging  sites  that  neighbor  the  same  surface  Al'^'^  (Ebamer  = 
27  kcal/mol).  However,  jumping  to  bridging  sites  on  a  different  surface  Al'^'^  atom  requires  a 
minimum  of  47  kcal/mol.  We  also  calculated  the  diffusion  barrier  for  a  concerted  diffusion  path 
in  which  displaces  an  which  then  moves  into  a  bridging  site  on  the  neighboring  Al'^'l 
The  barrier  to  concerted  O  diffusion  is  61  kcal/mol.  These  results  suggest  that  in  the  absence  of 
any  other  species,  O  becomes  localized  near  the  surface  Al'^'^  sites.  However,  several  of  the 
diffusion  paths  involve  singlet  to  triplet  conversion  and  need  to  be  further  studied  before  this 
conclusion  can  be  confirmed. 

Table  4.  Surface  diffusion  of  bound  O  and  H. 


Path  ID 

Oxygen 

Hydrogen 

H  assisted  O 

Ebarrier 

Ebarrier 

Ebarrier 

1 

4 

18 

- 

2 

4 

18 

- 

12 

22 

99 

- 

4 

27 

42 

1 

14 

27 

42 

3 

16 

47 

42 

3 

5 

46 

34 

45 

6 

47 

34 

46 

9 

67 

34 

43 

10 

65 

18 

48 

13 

64 

19 

48 

15 

64 

19 

46 

3 

68 

34 

- 

7 

49 

25 

- 

11 

68 

107 

- 

Note:  O  diffusion  along  paths  1,  2,  12,  3,  7,  and  1 1  involves  a 
triplet-singlet  spin  change  and  must  be  studied  in  more 
detail  before  the  diffusion  barriers  can  be  verified. 

We  also  studied  whether  the  presence  of  pre-adsorbed  H  can  promote  O  diffusion.  We  find  that 
pre-adsorbed  H  nearly  eliminated  the  barriers  to  O  diffusion  around  Al'^'^  sites  and  reduced  the 
barriers  along  paths  10,  13,  and  15  by  17  kcal/mol.  However,  this  improvement  was  not 
adequate  to  allow  O  to  completely  transverse  a  unit  cell.  Alternatively,  when  O  and  H  are 
allowed  to  diffuse  as  a  pair,  the  maximum  barrier  encountered  while  traversing  the  unit  cell  is  24 
kcal/mol.  These  results  suggest  that  the  presence  of  both  mobile  O  and  H  species  are  required 
for  O  surface  diffusion. 
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3.4  H2  Dissociation 


H2  cannot  molecularly  adsorb  to  the  A1  terminated  aAl203  surface;  however,  it  is  both 
thermodynamically  and  kinetically  favorable  for  H2  to  dissociate  into  1-4  (Ebarrier  =  20  kcal/mol, 
Erxn  =  -9  kcal/mol)  and  2-2  (Ebamer  =  2  kcal/mol,  Erxn  =  -60  kcal/mol)  dissociation  products. 
These  results  support  the  experimental  measurements  of  Wang  (1996),  who  postulated  that  H2 
adsorbed  to  the  AI2O3  support  resulted  in  a  dramatic  increase  in  the  H2  response  during  temporal 
analysis  of  products  (TAP). 

3.5  H  Surface  Diffusion 

From  the  results  in  table  4,  we  find  that  H  can  diffuse  from  a  Al'^'^  tetrahedron  site  (figure  4C)  to 
either  the  first  (Ebamer  =  1 8  kcal/mol)  or  second  nearest  neighboring  O  atom  (Ebarrier  =  26 
kcal/mol).  Once  in  these  conformations,  diffusion  between  sites  neighboring  the  same 
surface  Al'^*^  site  is  unlikely  (Ebarrier  =  42  kcal/mol).  However,  H  can  diffuse  between  sites 
neighboring  different  Al'^*^  sites  with  barriers  of  18  and  34  kcal/mol,  depending  on  how  close  the 
sub-surface  A1  atoms  are  to  the  diffusion  path. 


4.  Conclusions 


In  summary,  we  provide  a  theoretical  model  for  reaction  processes  that  occur  on  the  AI2O3 
substrate  and  influence  combustion  at  the  catalytic  particle.  Our  model  suggests  that  H2O  is  a 
primary  source  of  mobile  O  and  H  species.  Once  dissociated,  H  can  diffuse  away  from 
dissociated  H2O  leaving  behind  an  AI'^'Lqh  hydroxyl,  which  can  also  further  dissociate  creating 
a  2-2  0**H  pair.  Although  isolated  O  atoms  encounter  large  barriers  to  surface  diffusion,  the  2-2 
dissociated  0**H  pair  can  diffuse  with  a  barrier  ~24  kcal/mol.  We  also  find  that  H2  dissociation 
is  an  active  source  of  mobile  H  species.  Although  O2  can  adsorb  molecularly,  it  cannot  further 
dissociate  to  create  mobile  O.  However,  given  that  the  presence  of  pre-adsorbed  H  had  a 
profound  influence  on  the  surface  diffusion  of  O,  it  may  be  likely  that  H  can  also  influence  the 
dissociation  of  molecularly  adsorbed  O2. 
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